Frequency-dependent shot noise in long disordered SNS contacts 
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The shot noise in long diffusive SNS contacts is calculated using the semiclassical approach. 
At low frequencies and for purely elastic scattering, the voltage dependence of the noise is of the 
form Si = (4 A + 2eV)/?>R. The electron-electron scattering suppresses the noise at small voltages 
resulting in vanishing noise yet infinite dSi/dV at V = 0. The distribution function of electrons 
consists of a series of steps, and the frequency dependence of noise exhibits peculiarities at ui = neV, 
W = neV — 2A, and to = 2A — neV for integer n. 



PACS numbers: 72.70.+m, 74.40+k, 74.50+r 

In recent years much work was devoted to shot noise in 
hybrid superconductor (S) - normal-metal (N) systems.EJ 
The key role in these systems is played by Andreev reflec- 
tions at the NS boundaries, in which electrons incident 
from the normal metal are reflected as holes transferring 
a charge 2e to the superconductor. Most of theoretical 
results in this field were obtained under the assumption 
of quantum-coherent transport. The shot noise in diffu- 
sive SN contacts was shown to double at subgap voltages 
eV < A with respect to the noise in normal contacts,El 
which was interpreted as a doubling of the effective elec- 
tron charge. In SNS contacts, multiple Andreev reflec- 
tions from the contact edges-were found to greatly in- 
crease the low-voltage noiseETEl, which was attributed 
to independent transfers of large charge quanta 2A/eK 
times 2e. Both effects were observed experimentally.ETEl 
The frequency dependence of shot noise in SN contacts 
was shown to exhibit peculiarities at the Josephson fre- 
quency uj — 2eV at subgap voltageS|-and additional sin- 
gularities at uj — eV ± A at eV > An The peculiarity at 
u> = 2eV was already experimentally observed£3 

It was shown recently that the shot noise in diffusive 
SN contacts may be also described, using the semiclassical 
Boltzmann-Langevin approach.Ej In this approach, the 
doubled shot noise is explained by the effective heating pf 
the electron gas in the contact by the applied voltage.E3 
Because the SN boundary obstructs the heat flow from 
the contact into the superconducting electrode, the ef- 
fective heating is stronger and the noise is larger than in 
normal contacts. This approach also allows one to con- 
sider the effects of inelastic scattering on the noise. In 
this paper, the semiclassical approach is applied to the 
noise in long diffusive SNS contacts. 

Consider a narrow normal-metal microbridge connect- 
ing two massive superconducting electrodes. The elastic 
mean free path I of electrons in the microbridge is as- 
sumed to be much larger than (Z?/A) 1//2 , where A is the 
superconducting gap and D is the diffusion coefficient in 
the normal metal. We restrict ourselves to the case of 
zero temperature. The constant voltage bias is assumed 
to be much larger than the Thouless energy et = D/L 2 . 



Under these conditions, it is possible to neglect the pen- 
etration of the condensate into the microbridge and con- 
sider it just as a normal metal -with a nonequilibrium 
distribution function of electrons.EJ Since the SN bound- 
aries have a zero electrical resistivity and the whole volt- 
age drop takes place in the bulk of the contact, the 
noise in this systems is caused by the randomness of 
the impurity scattering and the semiclassical Bolizmann- 
Langevin approachMI or its quantum extensionlij may be 
applied. The Langevin equation for the current fluctua- 
tions reads 

Sj = -DWSp - <rVS(f> + Si ext , (1) 

where a is the electric conductivity, Sp(r) is the charge- 
density fluctuation, <5</>(r) is the local fluctuation of the 
electric potential, and the correlator of extraneous cur- 
rents Sj ext is given by 

Wr^O^O-a)}- = ^^(n - v 2 )T N (r 1 ), 

7V(r,w) = i J de{/(r,e)[l-/(r,e + u;)] 

+/(r,e + u;)[l-/(r,e)]}, (2) 

where /(r, e) is the average distribution function of elec- 
trons, which is almost isotropic in the momentum space 
because of a strong impurity scattering. If the frequency 
oj is smaller than 1 / RC, where R is the resistance of the 
contact and C is its capacity, the noise in the contact is 
obtained just by averaging Tn over its volume:E3 

Si=(i/R) J(dx/L)T N (x), (3) 

where x is the coordinate along the contact. To calculate 
the noise, one has to find first the average distribution 
function /, which obeys the standard diffusion equation 

DV 2 f + I in (e,x) =0, (4) 

with inelastic collision integral li n . What makes the sys- 
tem different from a normal contact is the boundary 
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FIG. 1. The scheme of energy bands for 2A/eV = 10/3. 
Dashed lines show the Fermi levels of the superconductors. 
The dashed areas show the path of a through passage of the 
quasiparticles across the contact. 

conditions for this equation. Consider first a single SN 
boundary located at x — 0. Assuming a total Andreev 
reflection at subgap energies, one may write the bound- 
ary conditions for the momentum-dependent distribution 
function at the SN interface in the formE3 



f(e,v x ) = 1 - f(-e,-v x ). 



(5) 



provided that the energy e is measured from the Fermi 
level of the superconductor. They are conveniently 
rewritten in terms of the symmetric and antisymmet- 
ric parts of f(s,v x ) with respect to v x . Expressing the 
antisymmetric part of / in terms of the gradient of the 
isotropic /(e, x), one obtains the following boundary con- 
ditions for Eq. (0) at the left end of the contact: 



/(e,0) = l-/(-e,0), 
V/(e,aO| a=0 = - V/(-e,aO| x=0 . 



(6) 



Making use of Eq. (|4|), one may express the gradient 
terms in (^|) in terms of the distribution function at the 
other end of the contact and the collision integral Ii ni 
which gives 



f(e, 0) = (1/2) [1 + /( e ,L) -/(- e ,L)] 



2D 



dx'(L - x')[I(e, x') - I(-e,x') 



(7) 



Now consider the case where both electrodes are super- 
conducting and their Fermi levels are shifted symmetri- 
cally by ieV/2. At first assume that the inelastic scat- 
tering is absent. Denote fie, 0) by and f(e,L), by 
/ij(e). Then one may write down the boundary condi- 
tions (Q) for the left and right interfaces in the form 



f L (e) = (1/2)[1 + f R (e) - f R (-e - eV)], 
f R (e) = (1/2)[1 + f L (e) - f L (-e + eV)] 



(8) 
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FIG. 2. The distribution functions at the right (solid line) 
and left (dashed line) ends of the contact for eV/2A = 3/10. 

Suppose first that the ratio 2A/eV — N is integer. In 
this case, the relevant range of energies —A — eV/2 < 
e < A + eV/2 may be split into N + 1 bands of width eV 
where the distribution functions //, and fn are constant. 
Denote the corresponding values by fn, and fm in order 
of increasing energy, i = 1 . . . N + 1. It follows from (^) 
and the continuity of / outside the subgap region that 
fiL and fm obey a set of equations 



fiL — (1 + fiR - f. 



N+l- 



i,fl)/2, 



i = 1 



.N 



fiR = (1 + fiL - /, 



N+3-i,L 



)/2, 



i = 2...N + l, 



/w+l,L — 0, 



fiR = 1. 



(9) 



It is noteworthy that at even N all equations (^) are 
coupled and a quasiparticle can traverse the contact via 
a series of Andreev reflections. For odd N this set of 
equations breaks into two independent subsets for odd 
and even i. In this case, the quasiparticles cannot tra- 
verse the contact via a series of Andreev reflections and 
are finally reflected back into the emitting electrode. 

Suppose now that the ratio 2A/eV is noninteger, i. e. 
2A/(N + 1) < eV < 2A/N. In this case, the relevant 
range of energies is split into alternating bands of width 
5i = 2A - NeV and S 2 = (N + l)eV -2A.N + 2 bands 
of width <5i and N + 1 bands of width 82 in all (see Fig. 
1). The values of the distribution functions in each set 

of bands />/£ RSj and f^l R ^ obey closed sets of equations 
(^) with numbers of bands N+2 and N+l. The solutions 
are of the form 

= 1 - i/(N + 2), /$ = 1 - (i - 1)/(JV + 2), 



f(2) _ 



1 - 



N + l 



f (2) = 1 _ 
JiR " 1 



i - 1 
N + l' 



(10) 
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A typical behavior of the distribution function for a non- 
integer ratio 2A/eV is shown in Fig. 2. It falls down 
from 1 to through a series of steps of unequal height. 
In principle, these steps ccwld be observed directly us- 
ing tunneling spectroscopy.£3 The distribution function 
at any point x is calculated from the diffusion equation 
as 

/(e, x) = (1 - x/L)f L {e) + (x/L)f R (e) (11) 

and integration in Eq. ([}]) is easily performed. At low 
frequencies, one obtains 

5/(0) = 2(2A + eV)/3R. (12) 

Despite the stepped distribution function, Si(V) is a per- 
fectly smooth curve. At large voltages, the SNS contact 
exhibits excess noise AA/3eV, whkjh is two times larger 
than that of a long SN contact .til At low voltages, it 
tends to the same finite value, as it takes-, place for a 
diffusive SNS contact with a strong barrier H The occur- 
rence of the same prefactor 1 /3 in the zero- voltage noise 
as in the voltage dependence of noise in normal-metal 
contacts is not purely coincidental. At low voltages the 
steps in the distribution function merge and it becomes 
a coordinate-independent linear function of energy in the 
range —A < e < A. On the contrary, the distribution 
function in normal-metal contacts is energy-independent 
linear function of coordinate in an energy range of width 
eV. Since the integration in the expression for the low- 
frequency noise is symmetric with respect to coordinate 
and energy, it gives the same prefactor in both cases. 

Though the stepped structure of the distribution func- 
tion does not manifest itself in the low-frequency noise, it 
results in peculiarities of its frequency dependence. The 
jumps of dSi I du) arise when the steps of /(e) come across 
the steps of f(e — u>) in the integrand of Eq. (||). There 
are three types of peculiarities of Si(lo). The first one 
corresponds to frequencies u> = meV, m — 1 . . . JV + 1. 
The magnitude of the jump is 

fl (m)_ 4 (JV + 2)(A+1)(2/V-3to + 3) + to 3 
3R (JV + l) 2 (JV + 2) 2 ' 1 j 

form= 1 . . . JV and 9 { Q N+1) = (<i/3R)(N+l)- 1 (N+2)- 2 . 
The second type corresponds to the frequencies ui — 
meV + 5% — 2A — (JV — m)eV, and the corresponding 
jumps of the derivative are 

(ro) _ 2 (JV + 2) (3mJV - 3m 2 - 2) + (JV + 2) 3 - to 3 
+ ~ 3R (JV + l) 2 (JV + 2) 

(14) 

for m = 1...JV and 6 { + +1) = (2/3R)(N + 2)~ 2 . The 
third type of peculiarities takes place at uj = meV — Si — 
(JV + m)eV — 2A resulting in the jumps of the derivative 

e {m) = ( 2 /3R)(N - m + 2) 3 /[(JV + 1)(JV + 2)] 2 (15) 
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FIG. 3. Frequency dependence of dSi/duj for eV/2A 
= 3/10 (solid line) and eV/2A = 1/3 (dashed line). 

for to = l...JVand ( ^ +1) = (2/3R)[(N+ 1)(JV + 2)]~ 2 . 
If 2A/eV is integer, Si = and the three sets of steps 
merge into one at u = meV . The magnitude of the steps 
is given by the sum of 9^ , 0+ , and 8^ 

= (4/i?)(JV+l-m)/(/V + l) 2 (16) 

for to = 1 . . .JV and 6 { * +1) = (2/3i?)(JV + l)" 2 . The fre- 
quency dependences of dSi/dto for an integer and nonin- 
teger ratios 2A/eV are shown in Fig. 3. 

Consider now the effects of electron-electron scatter- 
ing on the shot noise. This type of scattering slightly 
increases it in NN and SN contacts at eV A because 
it does not decrease the energy of electron-gas yet broad- 
ens the band of partially occupied states. M However the 
effects of electron-electron scattering in SNS contacts at 
eV <C A are essentially different from these cases since 
the overheated electron gas is trapped between two NS in- 
terfaces and quasiparticles can escape into the electrodes 
only by climbing to above-gap energies through a series 
of Andreev reflections or via inelastic scattering. Hence 
the electron-electron scattering should decrease the shot 
noise because electrons and holes outscattered from the 
subgap region immediately diffuse into the electrodes and 
the contact is thus effectively cooled. 

As the distribution function and the collision integral 
are only slightly coordinate-dependent at small voltages, 
one obtains from Eq. (f7j) 

f L (e) = (1/2) [1 + f R (e) - f R {-e - eV)} 



+ (L 2 /4£>) [7 in (e) - I in (-e - eV)}, 



f R (e) = (1/2) [1 + f L (e) - f L (-e + eV)} 
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FIG. 4. Plot of the dimensionless noise vs. dimensionless 
parameter of electron-electron scattering at eV <C A. 



one readily obtains that 5/ = AT e /R oc i? _1 A/ In A. This 
suggests that at finite a ee the noise tends to zero at zero 
voltage yet with an infinitely large slope. 

In summary, we have shown that in diffusive SNS con- 
tacts much longer than the coherence length the zero- 
temperature shot noise is a linear function of voltage with 
excess noise two times larger than in similar SN contacts. 
The finite zero- voltage noise 4A/3i? is suppressed by a 
finite electron-electron scattering, but the initial slope of 
Si(V) remains infinite. The stepped electron distribution 
function in the contact results in three series of peculiar- 
ities of the frequency dependence of noise at w = neV, 
lu = neV — 2A, and u> = 2A — neV for integer n. 



+ (L 2 /4D)[I m (e) - I m {-e + eV)]. 



(17) 



By expanding and fn to the second order and to 
zcroth order with respect to eV and forming symmen- 
tic and antisymmetric combinations of Eqs. ( p7| ) with 
repect to the indices I, R and ±e, one finally arrives at 
an energy-diffusion equation for the distribution function 
7(e) = [f L (e) + f R (e)]/2 in the form 

{eV) 2 d 2 J/de 2 + (L 2 /2D)[I(e) - I(-e)} = (18) 

with the boundary conditions /(—A) = 1 and /(A) = 0. 
The energy diffusion described by the right hand side of 
Eq. ([l8|) is the result of multiple Andreev reflections in 
the limit eV >C A. In the absence of inelastic scatter- 
ing it just gives f(e) = (1 — e/A)/2, which results in 
5/ = 4A/3i? in agreement with Eq. ([l2|). Making use of 
the explicit form of the collision integral 

I ee (s) = aeeep 1 J de' J duj{f(e - u)f{e') 
x[l-/(e)][l-/( £ '-u;)] 

-f(s)f(e' - W )[l - f(e - f(s')]} , (19) 

Eq. ( |i~8|) was solved numerically. The resulting dimen- 
sionless low-frequency noise RSi/AA is plotted in Fig. 
4 as a function of the dimensionless inelastic parameter 
A = a ee A 4 /(e 2 V 2 EpD/L 2 ). At large A the noise scales 
nearly as inverse logarithm of A. This behavior may be 
explained as follows. At large A the distribution function 
of electrons trapped in the contact assumes almost equi- 
librium shape with an effective temperature T e « A. On 
the other hand, the quasiparticles escape from the con- 
tact in this case mostly through electron-electron scat- 
tering rather than through Angreev reflections. Hence 
the Joule heating of the contact V 2 / R may be equated 
to the flux of energy J carried by electrons and holes 
outscattered from the subgap energy region by electron- 
electron collisions. As J oc a ee (A 3 T e /ep) exp(— A/T e ), 



1 Y. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 

2 M. J. M. de Jong and C. W. J. Beenakker, Phys. Rev. B 
49, 16070 (1994). 

3 D. Averin and H. Imam, Phys. Rev. Lett. 76, 3814 (1996). 

4 J. C. Cuevas, A. Martin-Rodero, and A. Levy Yeati, Phys. 
Rev. Lett. 82, 4086 (1999). 

5 E. V. Bezuglyi, E. N. Bratus, V. S. Shumeiko, and G. 
Vendin, Phys. Rev. Lett. 83, 2050 (1999). 

6 X. Jehl, P. Payet-Burin, C. Baraduc, R. Calemczuk, and 
M. Sanquer, Phys. Rev. Lett. 83, 1660 (1999). 

7 X. Jehl, M. Sanquer, R Calemczuk, and D. Mailly, Nature 
405, 50 (2000). 

8 T. Hoss, C. Strunk, T. Nussbaumer, R. Huber, U. Staufer, 
and C. Schonenberger, Phys. Rev. B 62, 4079 (2000). 

9 G. B. Lesovik T. Martin, and J Torres, Phys. Rev. B 60, 
11935 (1999); |cond-mat/0004489| (unpublished). 

10 A. A. Kozhevnikov, R. J. Shoelkopf and D. E. Prober, 
Phys. Rev. Lett. 84, 3398 (2000). 

11 K. E. Nagaev and M. Biittiker, |sond-mat/000712l| (unpub- 



lished) . 

12 By heating, we mean a broadening of electron distribution 
function, which does not necessarily imply thermalization 
of electrons. 

13 When the manuscript of this paper was in preparation, 
a preprint appeares where the low-frequency noise in this 
system was calculated, E. V. B ezuglyi, E. N. Bratus', 



:ond-mat/0010203 (un 



V. S. Shumeiko, and G. Wendin, 
published) . 

14 Yu. V. Nazarov and T. H. Stoof, Phys. Rev. Lett. 76, 823 
(1996). 

15 Sh. M. Kogan and A. Ya. Shul'man, Zh. Eksp. Teor. Fiz. 
56, 862 (1969) [Sov. Phys. JETP 29, 467 (1969)]. 

16 K. E. Nagaev, Phys. Rev. B 62, 5066 (2000). 

17 K. E. Nagaev, Phys. Rev. B 57, 4628 (1998). 

18 S. N. Artemenko and A. F. Volkov, Zh. Eksp. Teor. Fiz. 
72, 1018 (1977) [Sov. Phys. JETP 45, 533 (1977)]. 

19 H. Pothier, S. Gueron, N. O. Birge, D. Esteve, and M. 
Devoret, Phys. Rev. Lett. 79, 3490 (1997). 



4 



